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Foreword 


(U)  This  study  was  undertaken  to  examine  the  practicability  of  modeling 
the  effects  of  explosions  on  deep  underground  protective  structures,  with 
special  emphasis  given  to  modeling  the  response  of  unlined  cavities  in 
massive  rock  formations.  The  work  was  sponsored  by  the  Office,  Chief  of 
Engineers  through  the  U.  S.  Army  Research  and  Development  Project  "Military 
Engineering  Applications  of  Nuclear  Weapons  Effects  Research.  "  For  the 
most  part,  the  study  was  completed  during  the  summer  of  1964,  The  infor¬ 
mation  compiled  and  analyzed  in  this  report  was  used  as  a  basis  for  design¬ 
ing  a  small-scale  field  test  in  which  a  concrete  model  system  was  used  to 
simulate  the  effects  of  an  explosion  on  unlined  cavities  mined  within  a 
massive  homogeneous  rock,  (u) 

(U)  This  report  was  prepared  by  Dr.  George  B.  Clark.  The  author  desires 
to  acknowledge  the  helpful  suggestions  of  Mr.  John  N.  Strange  and  Lt.  A.  J. 
Hendron,  Jr.  ,  Chief  of  the  Engineering  Research  Branch  and  Military  As¬ 
sistant  to  the  Chief,  Nuclear  Weapons  Effects  Division  (NWED) ,  respectively. 
The  study  was  accomplished  under  the  overall  supervision  of  Messrs.  F.  R. 
Brown  and  G.  L.  Arbuthnot,  Chief  and  Acting  Chief  of  the  Nuclear  Weapons 
Effects  Division,  respectively,  (u) 

(U)  Col.  Alex  G.  Sutton,  Jr. ,  CE,  and  Col.  John  R.  Oswalt,  Jr. ,  CE. , 
were  .Directors  of  the  Waterways  Experiment  Station  during  the  preparation 
and  publication  of  this  report.  Mr.  J.  B.  Tiffany  was  Technical 
Director,  (u) 
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Glossary 

Sonic  velocity,  ft/sec 
Depth,  ft 

Closure  distance,  ft 

2 

Young's  modulus,  Ib/in. 

2 

Shear  modulus  of  elasticity,  lb/in. 

A  constant 
Length,  ft 
Mass,  lb 

Apparent  crater  radius,  ft 
Range,  ft 

Critical  range  at  which  crushing  ceases,  ft 
Pulse  duration,  msec 
Time,  sec 

Particle  velocity,  ft/sec 

3 

Apparent  crater  volume,  ft 
Yield  of  explosive,  lb,  kt,  or  Mt 
Depth  of  burst,  ft 
Strain,  dimensionless 
Radial  strain,  dimensionless 
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«u  Ultimate  strain,  dimensionless 

X  Model-to-prototype  length  ratio,  dimensionless;  also  Lame's 
(equation  2) 

X  Scaled  depth  of  burial  of  charge,  ft/lb or  ft/kt^’^ 

Xr  Scaled  range,  ft/lb1^  or  ft/kt1^ 
p  Model-to-prototype  mass  ratio,  dimensionless 
p  Mass  per  unit  volume  (mass  density),  lb/ft3 

p 

Rupture  stress,  lb/in. 

T  Model-to-prototype  time  ratio,  dimensionless 
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(u)  An  analysis  is  presented  of  the  scaling  parameters  which  are  impor¬ 
tant  in  scaling  explosive-induced  waves  in  earth  materials  and  the  effect 
of  these  waves  on  deep  underground  structures.  Basically,  two  approaches 
to  scaling  are  possible.  Gravitational  effects  can  be  allowed  for  and 
material  properties  scaled,  or  gravitational  effects  can  be  ignored  and 
material  properties  kept  the  same  in  the  model  as  in  the  prototype.  The 
study  indicates  that  peak  strain  is  dependent  upon  yield  as  well  as  other 
factors.  Closure  distances  resulting  from  confined  detonations  vary  from 
1,85  to  2.00  ft/lb1^,  the  variation  apparently  being  due  to  properties 
of  rock.  The  strain  magnitude  times  the  strain  pulse  period  together  with 
the  energy  level  of  the  pulse  appear  to  be  reliable  parameters  for  damage 
prediction.  High  explosives  (HE)  and  nuclear  explosives  (NE)  axe  believed 
to  be  almost  equal  in  effect  for  depths  of  burial  greater  than 

X  =0.20  ft/lb^^j  for  shallower  burial,  an  equivalence  factor  must 
c 

be  used.  (U) 

(u)  The  foregoing  analysis  indicates  that  model  scaling  can  be  used  to 
advantage  to  investigate  a  number  of  the  factors  which  are  important  in 
the  response  of  deep  underground  structures  to  explosive  attack,  (u) 
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SOME  BASIC  PRINCIPLES  OF  SCALING  EXPLOS ION- PRODUCED  DAMAGE 
TO  DEEP  UNLINED  OPENINGS  IN  ROCK  (U) 


Introduction 


Objectives 

(U)  1.  The  immediate  objective  of  the  study  reported  herein  was  to 
furnish  a  basic  foundation  for  investigating  the  feasibility  of  modeling 
some  of  the  important  parameters  affecting  damage  as  caused  by  surface 
or  near-surface  nuclear  detonations  to  deep,  unlined  protective  shelters 
in  rock.  The  design  of  model  experiments  and  instrumentation,  and  the 
development  of  standards  for  evaluating  the  test  results  will  require 
considerable  innovation.  However,  it  is  hoped  that  model  test  results  will 
furnish  a  basis  for  evaluating  quantitatively  the  effects  of  such  param¬ 
eters  as  geologic  structure,  cavity  size  and  shape,  and  various  liner  ma¬ 
terials  and  supports  on  the  overall  response  of  an  underground  inclusion 
to  explosions  in  the  megaton  range.  As  a  first  step  toward  developing  a 
modeling  scheme,  a  detailed  analysis  was  made  of  related  field  and  lab¬ 
oratory  experiments  and  the  data  obtained  from  them.  (U) 

Background 

(U)  2.  There  are  very  few  model  tests  of  record  which  deal  with  the 

1* 

destruction  of  underground  openings  by  explosives.  Jones  and  McCutchen 

conducted  tests  utilizing  black  powder  and  materials  with  properties 

scaled  according  to  the  laws  of  similitude ,  i. e.  a  gravity  ratio  of  1  was 

2  3 

assumed.  McCutchen  and  Arbuthnot  and  Strange  outlined  some  of  the  basic 
problems  involved  in  modeling  properties  of  materials  to  fit  selected  mod¬ 
eling  schemes.  (U) 

(U)  3,  Clark  and  Bruzewski,^  utilizing  natural  rock  materials  and  chem¬ 
ical  eaqplosives  in  tunnel  demolition  model  experiments,  showed  that  the 
cube-root  law  could  be  used  to  predict  closure  distances  for  charges  in 
the  range  of  1.0  to  6.0  lb.  Results  of  the  Underground  Explosion  Test 
(U£T)  program,  which  included  charges  weighing  up  to  320,000  lb  (TNT), 


*  Raised  numbers  refer  to  similarly  numbered  items  in  the  Selected 
Bibliography  at  the  end  of  this  report. 
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also  agreed  closely  with  cube-root  scaling  for  tunnel  closure  phenomena.  ** J ^ 

Several  very  small  scale  experiments  utilising  only  a  few  ounces  of  HE  have 
7 

been  reported,  but  the  results  for  the  most  part  are  too  inconclusive  to 
be  applicable.  In  addition  to  the  HE  tests  of  the  UET  program,  closure 
distances  are  now  available  for  several  underground,  contained  nuclear 
explosions  (table  l).  Limited  strain  and  particle-velocity  measurements 

Q 

were  obtained  from  a  few  of  the  experiments  listed  in  table  1.  (u) 

(C)  4.  In  general,  recent  analyses  of  cratering  information  indicate 
that  for  scaled  depths  of  burial  greater  than  approximately  0.20  ft/lb1^ 
(~25  ft/kt1^),  the  cratering  effectiveness  of  nuclear  explosives  (NE) 
and  high  explosives  (HE)  is  approximately  equal.  ^ 5 10  For  shallower  depths 
of  burial,  NE  become  increasingly  ineffective  compared  to  HE. 11  (c) 

Scaling  Theory  and  Analysis 
Approach  to  experimental  scaling 

(U)  5.  Any  method  of  experimental  scaling  of  explosion-produced  phe¬ 
nomena  depends  essentially  upon  the  application  of  appropriate  scaling 
laws  by  one  of  two  approaches:  (a)  by  utilizing  the  same  material  in  the 
model  and  the  prototype,  or  (b)  by  scaling  the  properties  of  the  prototype 
material.  In  most  explosion-effects  experiments,  as  a  matter  of  expedience 
the  same  materials  are  used  in  the  model  and  prototype.  If  nuclear  ex¬ 
plosions  are  to  be  modeled  at  small  scale,  it  is  necessary  to  model  the 
prototype  material  to  the  degree  practicable.  The  mechanical  effects  are 
much  the  same  for  HE  and  NE  at  confinements  exceeding  a  scaled  depth  of 
burial,  \q  ,  of  approximately  0.20  ft/lb1^.  For  confinements  less  than 
this,  a  calculated  equivalence  factor  relating  HE  and  NE  must  be  employed  to 
relate  model  HE  tests  to  full-scale  nuclear  explosions.  These  "adjusting" 
calculations  do  not  destroy  the  accuracy  of  experimental  results;  they 
do,  however,  affect  extrapolation  procedures.  Depending  on  the  modeling 
scheme  adopted,  there  are  other  parameters  that  may  require  empirical 
adjustment,  (u) 

(U)  6.  Generally,  the  effect  of  gravity  scaling  can  be  ignored  because, 
for  the  most  part,  gravity  has  no  direct  effect  on  underground  structures 
except  as  it  produces  lithostatic  pressures  which  can  be  artificially 
introduced  in  a  model.  In  a  given  experiment,  a  choice  must  be  made 
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between  (a)  allowing  for  gravitation  effects  and  scaling  material  prop¬ 
erties,  or  (b)  ignoring  gravity  and  keeping  material  properties  the  same 
in  the  model  as  in  the  prototype,  (u) 

Major  variables 

(U)  7.  A  list  of  parameters  which  might  well  govern  the  response  of  an 
underground  opening  or  inclusion,  along  with  the  scaling  relations  of  these 
parameters,  is  provided  in  table  2.  The  parameters  having  the  greatest 
influence  on  overall  response  must  be  selected  in  accordance  with  the  alms 
of  any  given  experiment.  The  paragraphs  that  follow  discuss  these  param¬ 
eters  and,  to  a  degree,  tell  how  they  relate  to  the  scaling  of  nuclear 
explosion  effects  on  deep,  underground  protective  structures.  (U) 

(U)  8.  Geometric  quantities.  The  quantities  angle,  area,  volume,  and 
curvature  are  independent  of  both  the  properties  of  the  material  and  grav¬ 
ity.  Thus,  if  all  lengths  are  properly  scaled,  geometric  scaling  is 
automatic.  (U) 

(U)  9.  Some  models  need  not  be  geometrically  similar  to  their  proto¬ 
types,  but  such  similarity  is  important  in  the  modeling  of  time-dependent 
phenomena  that  are  closely  related  to  overall  response.  In  most  cases 
time -dependent  phenomena  are  important  parameters  of  explosion-effects 
modeling.  Also,  it  is  not  always  practicable  to  model  all  features  of 
a  prototype  in  absolute  detail,  and  it  is  generally  accepted  that  only 
features  that  affect  the  response  of  a  structure  to  a  measurable  degree 
need  be  modeled.  For  example,  it  is  not  necessary  to  model  every  struc¬ 
tural  detail  of  an  air-vent  opening's  connection  to  a  deep  underground 
structure  in  order  to  determine  the  collapse  pressure  of  the  structures 
main  section,  (u) 

(U)  10.  In  modeling  a  large  rock  mass,  it  is  necessary  for  economic 
reasons  to  limit  the  lateral  dimensions  of  the  model;  whereas  a  nuclear 
event  of  interest  may  occtir  at  or  near  the  surface  of  a  solid  half-space. 
One  primary  concern  in  the  design  of  the  model  is  to  make  it  large  enough 
to  record  the  effects  of  the  incident  pulse  before  the  effects  of  re¬ 
flected  waves  are  registered  (superimposed)  upon  the  incident  waves'  sig¬ 
nature  or  response.  (U) 

(U)  11.  In  homogeneous  elastic  material,  the  stress  concentration  about 
a  circular  opening  is  independent  of  the  opening's  size;  however,  the 
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ratio  of  tunnel  diameter  to  the  effective  length  of  the  strain  pulse  deter¬ 
mines  the  response  mechanism  and,  very  likely,  the  extent  of  failure.  (U) 

(U)  12.  Frequency.  Frequency  has  the  dimensions  of  inverse  time  and 
should  scale  as  the  inverse  of  scaled  time.  This  means  that  the  rise  and 
decay  times  will  likewise  decrease  with  the  model  size  and  will  thus  have 
an  effect  upon  the  mechanism  of  failure  within  closure  distances  (see  sec¬ 
tion  on  "Strain  and  wavelength,"  paragraphs  31-41 ).  When  wavelengths  are 
large  compared  to  the  opening  dimension,  fracture  will  occur  by  crushing; 
when  relatively  short  wavelengths  prevail  as  compared  to  the  dimension  of 
the  opening,  the  response  will  likely  take  the  form  of  tensile  slabbing.  (U) 
(U)  13.  Velocity.  If  velocity  attenuation  due  to  material  properties 
is  ignored,  all  velocity  parameters  (shock,  sonic,  and  particle  velocity) 
are  the  same  in  the  model  and  prototype  when  the  model  and  prototype  mate¬ 
rials  are  the  same.  This  directly  implies  that  the  strain  will  also  be 
identical,  as  indicated  by  the  equation 

v  =  ce  (l) 

where  v  is  particle  velocity,  c  is  sonic  velocity,  and  e  is  the 
strain.  (In  equation  1,  if  there  is  a  difference  between  the  velocity  of 
the  peak  stress  or  strain  rates  of  propagation  and  the  sonic  velocity, 
velocity  of  the  peak  stress  should  be  used  fo*’  c  and  not  the  sonic 
velocity.  )  Sonic  velocity  is  given  by  the  equation: 


where  c  is  sonic  velocity,  X  is  Lame’s  constant,  G  is  the  shear  mod¬ 
ulus  of  elasticity,  and  p  is  the  mass  density.  Therefore,  equation  1 
will  hold  when  material  properties  are  modeled,  (U) 

(u)  14.  Accelerations.  If  gravitational  acceleration  is  ignored  (it 
should  not  be  when  the  gravity- induced  stresses  are  a  significant  part  of 
the  medium's  strength),  the  accelerations  associated  with  the  stress  wave 
will  then  scale  inversely  as  the  model  ratio  if  inputs  are  properly  scaled, 
i. e.  the  smaller  the  model,  the  higher  the  acceleration.  Ignoring  gravita¬ 
tional  acceleration  is  normally  Justified  because  gravity  affects  only  two 
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parameters  related  to  nuclear  attack.  The  first  is  the  formation  of  the 
apparent  crater  at  the  surface,  which  is  only  indirectly  related  to  the 
transmitted  energy.  The  second  is  lithostatic  pressure,  which  is  a  crit¬ 
ical  factor  in  determining  the  necessary  depths  for  reducing  damage  to 

on  acceptable  level.  This  may  he  demonstrated  as  follows.  Utilizing 

12 

Newmark's  eq^tion  for  stress  as  related  to  yield  and  range, 

Or,25(w)5/6(^0  )5/2(_c_)2  (3 

where 

r>. 

a  =  rupture  stress,  lb/in. e‘ 

W  =  yield  of  explosive,  Mt 
R  =  range,  ft 
c  =  sonic  velocity,  ft/sec 

p 

and  assuming  that  the  lithostatic  pressure  is  1  Ib/in.  per  foot  of  depth, 
the  combined  effects  of  dynamic  and  static  stress  for  three  yields  indi- 

P 

cate  that  for  a  rock  of  25,000  lb/in.  strength  there  is  no  depth  at 
which  unlined  tunnels  would  survive  in  a  100 -Mt  attack  (fig.  l).  Equa¬ 
tion  3  has  been  modified  to  fit  the  results  from  Project  HAKDHAT  in  combi¬ 
nation  with  an  analysis  of  UET  data  (see  "Strain  and  wavelength,"  para¬ 
graphs  31-41).  (U) 

(U)  15.  Generally,  models  do  not  reproduce  the  scaled  effects  of  litho¬ 
static  pressures  and  are  therefore  subjected  to  dynamic  stresses  only. 

It  is  possible,  however,  to  simulate  the  effects  of  lithostatic  pressures 
in  a  model  by  using  a  suitable  framework  and  heavy  springs  for  loading  to 
the  desired  static  stress  level;  or,  for  a  small  model,  a  centrifuge  may 
be  used  to  properly  simulate  these  stresses.  (U) 

(U)  16.  Angular  velocity  and  angular  acceleration.  These  quantities 
are  important  parameters  only  when  revolving  bodies  are  being  studied; 
thus,  they  have  no  bearing  on  the  problem  being  considered  herein.  (U) 

(U)  17 «  Density.  When  the  same  materials  are  utilized  in  model  and 
prototype,  the  density  ratio  is  specified,  and  when  one  other  ratio  is 
fixed,  all  other  model  ratios  are  determined.  When  different  materials 
are  used  in  model  and  prototype,  a  relatively  wide  choice  of  density 
ratios  exists.  Since  compressibility  of  the  materials  is  often  important 
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Fig.  1. 
dynamic 


The  effect  of  superimposing  gravity- induced  stresses  o 
free-field  stress  curves  for  surface  bursts  from  weapon 
yields  of  1,  10,  and  100  Mt 
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in  determining  the  net  response  of  a  system,  it  is  quite  unlikely  that  the 
density  ratio  will  exceed  5  (prototype  density  divided  by  model  density),  (u) 
(U)  18.  Momentum,  moment  of  momentum,  and  angular  momentum.  Of  these 
three,  only  the  momentum  is  considered  to  be  significant  in  the  problem 
under  study.  An  important  consideration  in  evaluating  damage  and  response 
effects  is  the  conservation  of  momentum.  It  is  scaled  as  shown  in  table  2.  (U) 
(U)  19.  Force.  Force  scales  as  the  mass  or  as  the  model  ratio  squared, 
depending  on  whether  or  not  the  materials  are  scaled.  However,  stress  and 
pressure  are  the  most  significant  force  parameters  of  concern  in  this  study,  (u) 
(U)  20.  Total  impulse.  This  parameter  varies  as  shown  in  table  2.  While 
it  has  not  been  demonstrated,  there  is  some  indication  that  the  peak  stress 
(or  strain)  is  not  the  most  important  factor  in  producing  damage.  The 
duration  and  magnitude  of  peak  stress  (or  strain)  above  a  certain  critical 
level,  i.e. ,  the  equivalent  of  total  impulse  above  the  critical  level  of 
Stress  or  strain,  appear  to  be  the  most  important  damage  parameter.  Im¬ 
pulse  follows  the  cube-root  law  of  scaling  for  similar  materials,  and  may 
account  in  part  for  the  fact  that  closure  distances  also  follow  this  law 
(see  "Strain  and  wavelength,"  paragraphs  31-^1).  (u) 

(u)  21.  Torque,  power,  and  action.  These  quantities  scale  as  shown  in 
table  2,  but  have  no  direct  application  to  the  general  problem  being  con¬ 
sidered  in  this  study. 

(C)  22.  Work  and  energy.  These  quantities  vary  as  the  cube  of  the  model 
ratio  for  similar  materials  and,  hence,  follow  the  cube-root  law.  When 
materials  are  modeled  (gravity  and  acceleration  ratios  equal  l),  the  fourth- 
root  law  is  applicable,  (c) 

(U)  23.  In  the  modeling  of  a  large  nuclear  explosion,  two  problems  arise: 

(a)  determination  of  the  dissipation  of  energy  in  the  spherically  expanding 
wave,  and  (b)  determination  of  the  equl valence  between  effective  energy 
yield  from  NE  and  HE,  particularly  for  shallow  depths  of  burst.  Inasmuch 
as  model  studies  must  employ  chemical  (HE)  explosives,  the  similarities 
and  differences  between  the  two  types  of  explosives  must  be  evaluated. 

Early  studies  of  the  cratering  effects  of  nuclear  explosives  led  to  the 
erroneous  conclusion  that  chemical  explosives  are  always  much  more  effi¬ 
cient  than  nuclear  explosives.  Analysis  of  more  recently  acquired  crater¬ 
ing  data  indicates  that  both  types  of  explosive  have  approximately  the 
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same  crater-producing  capability  when  scaled  depths  of  burial,  >v  ,  are 

greater  than  0.20  ft/lb1'3  (~25  :  .  t3^3).  Vaile11  utilized  curves  which 

assumed  equal  efficiency  below  a  A«pth  closely  approximating  that  given 

above  (X  =  0.14  ft/lb 3-'3) ;  and  Vortman^  showed  that  both  are  approximately 
c 

equal  in  producing  craters  for  conditions  of  deep  burial.  Other  calcula- 

13 

tions  which  resulted  in  low  NE  efficiencies  appear  to  be  based  on  the 
false  assumption  that  there  are  shock  losses  in  NE  but  not  in  HE,  and  that 
the  efficiency  of  the  latter  is  dependent  only  on  coupling.  The  use  of 
1/3.4  scaling  or,  in  fact,  any  type  of  scaling  which  will  normalize  crater 
parameters  to  a  common  curve  signifies  that  all  charges  of  that  size, 
whether  HE  or  HE,  are  of  equal  efficiency  with  respect  to  the  phenomena 
considered  and  the  given  parameter,  (u) 

(C)  24.  For  reliable  model  investigations,  it  is  particularly  necessary 

to  know  the  relative  efficiencies  of  shallow  HE  and  NE  as  compared  to 

completely  confined  detonations.  The  methods  used  by  Vailex  offer  few 

possibilities  because  of  the  uncertainty  of  choosing  a  specific  and  proper 

scaling  law  that  is  applicable  to  shallow  explosions.  (C) 

(C)  25.  If  it  is  assumed  that  for  scaled  depths  (cube-root  scaling)  of 

burial  less  than  Xc  =  1.0  ft/lb1/3  the  apparent  (or  true)  crater  volume, 

radius,  or  depth  is  a  measure  of  the  shock  energy  input  into  the  earth 

below  the  crater,  then  approximate  relative  efficiencies  can  be  determined. 

.he  energy  directed  downward  into  the  earth  approaches  a  maximum  at  a 

1/3 

burial  depth  of  X  =3.0  ft/lb  '  D  ,  and  cratering  parameters  approach  a 

^  l/o 

maximum  at  1.0  ft/lb  ' 0  depending  upon  charge  size.  That  is,  when  cube- 
.*oot  scaling  is  used,  the  maximum  values  of  the  scaled  radius,  volume,  and 
depth  tend  to  decrease  with  increase  in  yield.  (C) 

l)  26.  Apparent  crater  radii  in  desert  alluvium,  which  approximate  the 
ue  crater  radii  in  this  medium,  vary  with  yield  and  scaled  depth  as 
shown  in  fig.  2.  The  curves,  based  upon  cube-root  scaling,  are  shown  as 
two  straight-line  segments.  These  curves  were  derived  from  the  data 
given  in  table  3  and  are  for  radii  at  equal  scaled  depths  of  burial  over 
the  HE  ranges  from  256  to  40,000  lb,  and  for  NE  yields  up  to  100  kt.  The 
position  of  the  Xc  =  0  (HE)  line  was  determined  largely  on  the  basis 
of  the  JAHGLE  S  and  J0HNIE  BOY  events.  The  relative  radius -producing 
efficiency,  easily  determined  for  any  two  events  at  a  scaled  depth  of 
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Fig.  2.  Variation  of  scaled  apparent  crater  radius  with 

of  weapon  yield 
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1.0  ft/lbx'J,  can  also  be  considered  a  measure  of  the  relative  strain- 

producing  efficiency.  These  values  for  various  scaled  depths  of  burst 

were  utilized  in  connection  with  the  strain -producing  efficiency  for  HE 
14 

in  granite  as  a  basis  for  predicting  the  effects  of  very  large  HE  or 
NE  charges  (fig.  3).  (C) 

(C)  27.  The  data  in  table  4  (all  from  events  having  yields  of  from  0.5  to 
1.2  kt  except  SEDAN  which  had  a  yield  of  100  kt)  show  that  radii  efficien¬ 
cies  are  not  a  good  measure  of  relative  efficiencies  for  above-surface 
shots  (JANGLE  S),  but  are  appropriate  for  depths  of  burial  as  deep  as 

X  =  1  .  The  plot  of  scaled  crater  radius  versus  yield  (fig.  2)  indicates 
c 

that  the  crater  radius-producing  ability  of  NE  decreases  with  increasing 

weapon  yield  approximately  as  that  of  HE  except  at  shallow  burials,  where 

the  specific  crater  radius  is  smaller  for  NE.  This  percentage  of  energy 

12 

efficiency  is  somewhat  larger  than  that  suggested  by  Newmark.  The  crater 
radius -producing  efficiency  of  HE  is  much  higher  than  that  of  NE  for  deto¬ 
nations  both  immediately  above  and  below  the  surface,  (c) 

(C)  28.  For  burial  depths  greater  than  X  =0.20  ft/lb1^  ,  the  ef- 

c 

ficiency  of  NE  rapidly  approaches  that  of  HE.  This  is  supported  by  the 
following  conclusions  from  references  9  and  10: 

a.  The  apparent  crater  parameters  for  both  NE  and  HE  detona¬ 
tions  are  approximately  the  same  when  the  scaled  depths  of 
burial  are  equal  (fig.  2).^ 

b.  Peak  strains  from  completely  confined  NE  and  HE  explosions 
of  equal  yield  follow  approximately  the  same  law  at  equal 
scaled  distances. 

c.  An  analysis  of  the  gross  thermodynamic  processes  taking 

place  within  the  final  cavity  limits  produced  by  both  types 

of  explosive  indicates  that  the  heat  losses  in  NE  processes 

may  be  less  than  10  percent  greater  than  those  in  HE,  the 

former  being  due  to  the  heats  of  fusion  and  vaporization 

of  rock  etc.  ,  which  represent  relatively  small  irrecoverable 
10 

energy  losses. 

d.  Shock  losses  associated  with  HE  and  NE  should  be  approximately 
equal. 

e.  The  cube  root- scaled  closure  distance  for  HAEDHAT  (an  NE 
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SCALED  DEPTH  OF  BURIAL,  FT/LB1'* 


Fig.  3*  Percent  strain  produced  relative  to  completely  contained 
charges,  based  on  extrapolation  of  curve  from  reference  15  to¬ 
gether  with  crater  radius  curves  versus  yield  for  desert  alluvium 
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experiment)  was  greater  than  for  UET  KB  experiments  even 
allowing  for  the  relatively  shallow  burial  (Xc  a  0.36  lb/ft  ^^) 
of  the  UET  detonations. 

f.  With  the  exception  of  the  strain  and  velocity  pulse  lengths 
from  SKOAL,  the  NE  and  HE  values  obey  a  common  scaling  law 
for  UET  and  HAEDHAT  when  adjustments  are  made  for  the  charge 
size. 

£.  Subsidence  craters  created  by  completely  confined  NE  and  HE 
obey  a  common  cube-root  law,^  (C) 

(C)  29.  Thus,  for  experimentation  with  models,  chemical  explosives  may 
be  used  to  simulate  nuclear  explosives  on  a  1:1  basis  for  confinements 
varying  from  complete  to  scaled  depths  of  burst  approaching  a  minimum  value 
of  Xc  =  0. 20  ft/lb"^  .  For  shallower  depths  an  equivalence  factor  must 
be  used.  (C) 

(C)  30.  Stress,  pressure,  and  strength.  Theoretically,  these  quantities 
all  have  1:1  values  at  the  same  scaled  range  in  a  model  composed  of  the 
same  material  as  the  prototype.  However,  if  the  prototype  is  scaled  to 
the  megaton  range,  the  stress  at  given  scaled  ranges  from  the  explosive 
will  be  decreased  by  attenuation  factors  not  otherwise  accounted  for  by 
similitude.  Thus,  model  results  should  give  higher  values  for  dynamic 
stress  than  the  prototype;  however,  this  may  be  more  than  compensated  for 
by  the  stress  due  to  lithostatic  or  tectonic  pressure  (see  figs.  4  and  5). 

The  attenuation  oan  be  accounted  for  by  utilizing  an  expression  which  in¬ 
cludes  attenuation  (see  "Strain  and  wavelength"  below).  When  properties 
of  materials  are  modeled  but  with  density  ratio  equal  to  1,  then  the 
stresses  and  strengths  vary  as  the  first  power  of  the  model  ratio.  (C) 

(C)  31.  Strain  and  wavelength.  The  peak  strain  has  been  justifiably 
regarded  by  Newmark^  as  one  of  the  more  important  design  parameters  for 
tunnel  linings  and  packing  behind  the  liners.  The  crushing  of  rock  around 
a  tunnel  liner  is  assumed  to  be  accompanied  by  bulking,  and  the  range, 

,  at  which  crushing  in  granite  will  cease  is  given  for  a  contained 
burst  by  the  equation 

r =  °-°°33  =  °-01  ln-/ln-  ( rm )5/6  ( rf  )5/a  w 
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LITHOSTATIC  PRESSURE 
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LITHOSTATIC  PRESSURE 

EFFECTIVE  FAILURE  STRAIN  WITH 

LITHOSTATIC  PRESSURE 

EFFECTIVE  FAILURE  STRAIN  WITHOUT 

LITHOSTATIC  PRESSURE 

EFFECTIVE  FAILURE  IMPULSE  WITHOUT 

LITHOSTATIC  PRESSURE 

ADDITIONAL  EFFECTIVE  IMPULSE  WITH 

LITHOSTATIC  PRESSURE 


TIME 


g.  b.  Schematic  of  stress  or  strain  versus  time  illustrating  the  effects 

of  lithostatic  pressure 
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Fig.  5*  Schematic  showing  how  total  impulse  is  maintained  with  increase 
in  pulse  length  and  decrease  in  peak  dynamic  strain,  the  latter  being 
increased  by  the  lithostatic  pressure  with  increased  weapon  yield  and 

depth  of  underground  structure 
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where  eu  is  ultimate  strain W  is  the  explosion  yield  in  Mt,  and 
Rc£  is  the  critical  range  at  which  crushing  ceases.  Beyond  the  distance 
Rc(j  ,  there  will  be  some  crushing  and  failure  of  zone  2  severity  (contin¬ 
uous  spalling  but  with  decreasing  thickness),  (c) 

(C)  32.  The  net  response  is  thus  related  to  the  length  of  the  strain 
pulse,  which  has  been  shown  to  increase  approximately  as  the  first  power 
of  travel  distance.  Also,  the  pulse  duration  is  a  function  of  yield  and 
is  given  by  the  following  empirical  equation  which  was  derived  from  fig.  6. 

t  =  0.038  W0,4  ^  -5^  J  (5) 

where  t  is  pulse  duration,  R  is  range,  and  W  is  the  explosion  yield  in 
pounds.  (C) 

(C)  33.  The  only  available  strain  data  for  very  small  charges  are  for 

17 

Kanawa  sandstone  and  Lithonia  granite  '  which  have  properties  similar  to 
Navajo  sandstone  and  HARDHAT  granite.  However,  charges  were  fired  at  var¬ 
ious  scaled  depths,  and  hence  the  curve  for  S1.—II  charges  is  included  for 
reference  only,  (c) 

(C)  34.  Another  condition  which  increases  the  effective  strain  is  the 
lithostatic  pressure,  i.e.  ,  the  rock  is  already  under  vertical  compressive 
strain  due  to  the  weight  of  the  overlying  rock.  This  is  illustrated  in 
figs.  1,  4,  and  5  which  indicate  that  the  lithostatic  pressure  not  only 
increases  the  magnitude,  but  also  the  duration  of  the  effective  strain,  (c) 
(C)  35.  Thus,  the  values  of  strain  and  pulse  duration  should  be  increased 
by  a  factor  which  is  a  linear  function  of  depth  ,  d  : 


In  equation  6,  e  is  the  strain  (dimensionless),  k  is  a  constant,  W  is 
the  explosion  yield  in  Mt,  R  is  the  range  in  feet,  and  d  is  the  depth 
for  a  given  observation  in  feet.  In  equation  7>  t  is  the  shock  duration 
in  msec,  k  is  a  constant,  W  is  the  weight  of  charge  in  pounds,  and  R  is 
the  range  in  feet.  These  depth  factors  may  be  at  least  partially  responsible 
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for  the  fact  that  the  HARDHAT  closure  distance  was  larger  than  predicted. 

A  further  modification  of  the  above  strain  equation  is  developed  in  para¬ 
graph  3?.  (C) 

(C)  36.  The  maximum  strain  that  the  rock  medium  will  sustain  without 
failure  and  thus  transmit  to  the  rock  located  in  the  direction  of  propaga¬ 
tion  is  determined  by  the  breaking  strain  of  the  rock  under  the  particular 
conditions  of  confinement  at  the  time  it  is  traversed  by  the  pressure  wave. 
This  strength,  in  turn,  also  defines  the  limit  of  the  fracture  zone.  The 
properties  of  rock  in  the  UET  and  HARDHAT  programs  are  listed  in  table  5.  (C) 

(C)  37.  Still  other  factors  which  may  partially  account  for  the  greater 
scaled  closure  distance  for  HARDHAT  are  the  lower  strain  at  rupture,  and 
the  fact  that  there  is  less  than  half  the  energy  under  the  stress-strain 
curve  at  failure  for  HARDHAT  granite  as  for  UET  (Unaweep.)  granite.  If  the 
rock  is  elastic,  the  higher  Young's  modulus  for  HARDHAT  granite  would  have 
the  effect  of  increasing  the  seismic  velocity,  thus  decreasing  the  length 
of  the  pulse.  Generally,  the  slope  of  the  peak  strain  versus  scaled  range 
curve  for  granite  is  taken  to  be  -2. 5.  The  plot  of  points  for  320-lb  shots 
at  X  =  3*64  ft/lb  ,  which  is  assumed  to  approach  complete  confinement, 
and  the  peak  strains  for  HARDHAT  and  SHOAL  granite  with  RANTER  tuff  data 
added  for  comparison  (fig.  7) ,  and  similar  curves  for  sandstone  (fig.  8), 
give  the  following  peak  strain  equation: 

e  =  0.08  W"1/9  ^  J  ^  (8) 


where  the  constant  0. 08  has  dimensions  of  Vr'*  ,  W  is  in  pounds,  and 
R  is  in  feet.  That  is,  the  peak  strain  is  a  function  of  charge  size  as 


well  as  seeled  range.  The  data  from  SHOAL  were  not  employed  to  obtain  the 
multiplier  for  the  r/w1/3  term,  and  the  exponent  of  W  would  be  a 
little  smaller  than  1/9  if  these  were  included.  Ihe  strain  equation 
thus  becomes  an  expression  for  a  family  of  curves,  and  this  fact  should 
be  kept  in  mind  when  the  equation  is  employed  for  prediction  purposes. 
Equation  8  for  conf ined  explosions  may  be  converted  to  the  form  used  by 
Heumark  for  range  in  thousands  of  feet  and  yield  in  megatons: 
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R/WV3  -  SCALED  RANGE,  FT/LB1/3 

Fig.  8.  Variation  of  peak  strain  in  Navajo  sandstone  with 
scaled  range  for  specific  charge  weights  (Xc  =  0.36) 
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(C)  38.  It  is  also  assumed  that  the  energy  effects  of  NE  and  HE  are 

equal  for  confined  detonations  and  for  buried  detonations  below  a  scaled 

depth  of  X  =0.20  ft/lb  ^/3  .  For  near- surface  events,  the  equivalence 
c 

can  be  closely  approximated  as  follows:  The  crater  radius  appears  to  be 
a  fair  measure  of  the  relative  amount  of  energy  transmitted  to  the  earth. 
The  plot  of  radii  for  significant  HE  and  NE  events  (fig.  2)  also  indicates 
that  for  contact  bursts,  NE  is  approximately  60  percent  as  effective  as 
HE.  Utilizing  also  the  effect  of  depth  as  measured  from  the  UET  320-lb 
shots,  the  strain  energy  from  near-surface  NE  explosions  can  be  approxi¬ 
mated  (fig.  3)-  It  is  found  that  the  equivalence  for  1  kt  at  the  surface 
is  about  9*5  percent,  and  that  for  1  Mt  about  6  percent,  the  percentage 
decreasing  by  about  1  percent  for  each  order  of  magnitude  increase  in 
yield.  aus  for  the  1-Mt  range,  the  strain  equation  for  a  surface  burst 
becomes 

e  =  0.004(0.06  W)2/3  ^  ^  (10) 

or 

€ = 0. 0006(H)2/3  (n) 

12 

This  compares  f,  'ably  with  the  equation  recommended  by  Newmark  with 
the  exception  e-  the  exponent  of  the  yield  which  is  2/3  instead  of  5/6 

(fig.  9).  (c) 

(C)  39.  Young's  modulus  (E)  for  HARDHAT  granite  is  11.8  X  10^  lb/in. 2 
and  the  seismic  velocity  equals  18,000  ft/sec.  The  stress  equation  thus 
becomes 

CTr  -  Vs  -  36-400  x  o-oc06  w273  (  ^i22  )5/2  ( 5300  y 

(12) 

ar  -  22  W273  (  3§22)  5/2  (jgjjj)2 
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Fig.  9.  Effect  of  superimposing  gravity-induced  stresses  on 
dynamic  free-field  stress  curves  for  weapon  yields  of  1,  10, 
and  100  Mt  (surface  bursts).  Differs  from  fig.  1  because  of 
the  use  of  the  modified  stress  equation  12 
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Again,  W  is  in  Mt,  R  is  in  feet,  and  c  is  in  ft/sec;  is  radial 

strain.  The  stress  equation  is  also  logically  in  close  agreement  with 

Newmark  with  the  exception  of  the  yield  exponent.  (C) 

(C)  4o.  The  reduction  of  the  yield  exponent  from  5/6  to  2/3  reduces 

predicted  peak  stress  for  multimegaton  bursts.  Newmark^  maintained  that 

the  peak  strain  is  the  most  reliable  criterion  for  design  purposes. 
l8 

Clark-1  proposed  that  both  the  magnitude  and  the  duration  of  the  pulse 
above  a  critical  level  are  important  along  with  the  related  energy  density. 

A  plot  of  all  available  data  for  closure  distances  versus  yield  shows  that 
values  all  fall  within  a  scaled  range  of  1. 29  to  2.00  (fig.  10).  Excluding 
tuff,  the  range  is  1.65  to  2.00.  If  UET  results  are  adjusted  to  account 
for  the  fact  that  burial  depths  were  only  0.  36  ft /lb1/3 ,  and  ERDL  results 
in  basalt  are  adjusted  to  account  for  the  low  loading  density,  then  all 
closure  distances  (except  in  tuff)  for  large  yields  (>750  lb)  fall  very 
close  to  2.00  ft/lb1/3.  Thus,  the  scaled  closure  distance,  d  ,  in  feet 
(fig.  10)  can  be  expressed  mathematically  as 

dc  =  k  W1/3  (13) 

If  the  peak  strain  is  multiplied  by  the  strain  period  the  following  results: 


e  x  t  =  kQ  W1/3  (Xr)“1,5 

That  is,  at  a  given  scaled  range,  X^  ,  the  product  of  the  peak  strain  and 
pulse  duration  varies  approximately  as  the  cube  root  of  the  yield.  This 
is  analogous  in  form  with  the  closure  equation  (equation  13),  which  indi¬ 
cates  that  a  more  basically  correct  design  criterion  is  the  product  of 
the  strain  and  its  duration.  The  later  is  assumed  to  be  proportional  to 
the  area  under  bhe  strain-time  curve.  (C) 

(C)  4l.  There  is  no  immediately  obvious  explanation  for  the  smaller 
closure  distances  occurring  in  tuff.  The  lithostatic  pressure  would  be 
almost  comparable  to  that  of  HARDHAT.  Hence,  the  greater  resistance  of 
tuff  to  closure  must  be  attributed  to  its  physical  properties.  Two 
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properties  th&t  are  radically  different  from  those  of  the  other  rocks 
studied  are  the  compressive  strength  and  Young's  modulus.  Thus,  while  the 
magnitude  of  strain  In  the  RAINIER  wave  was  larger  than  that  of  HARDHAT, 
the  relative  "elastic"  energy  would  have  been  l/2  X  1500  x  0. 007  =  0.525 
for  tuff  compared  to  1/2  X  1500  x  0.002  «  12.5  for  HARDHAT  granite,  or 
approximately  25  times  as  great  for  the  latter.  While  these  figures  are 
for  unconfined  rock,  they  are  strongly  indicative  of  the  reason  for  the 
resistance  of  tuff,  or  conversely,  its  inability  to  sustain  waves  at  high 
energy  and  strain  levels.  (C) 

(U)  42.  In  a  model  wherein  the  properties  of  the  material  are  scaled, 

l/2 

the  impedance  varies  as  pX  '  where  X  is  the  ratio  between  prototype 

and  model  length;  and  if  the  densities  are  the  same,  the  impedance  will 
l/2 

vary  as  X  '  or  as  the  velocity  ratio.  In  a  model  which  has  the  same 
properties  as  the  prototype,  the  impedances  of  the  two  are  identical.  In 
either  case,  impedance  matching  requires  that  pc  of  both  model  and 
prototype  explosions  and  materials  should  be  scaled  according  to  the  se¬ 
lected  model  criteria,  (u) 

(C)  43.  In  the  calculation  of  impedance  matching  of  explosives  to 

15  19-22 

rock,  the  Bureau  of  Mines  uses  detonation  velocity  as  the  imped¬ 

ance  velocity  parameter.  Inasmuch  as  the  detonation  state  exists  only 
for  a  very  small  period  of  time  relative  to  the  effective  pressure  pulse 
of  a  confined  explosive,  it  would  appear  more  logical  to  employ  the  density 
and  sonic  velocity  of  the  explosion  state.  A  comparable  parameter  for 
nuclear  explosives  is  not  so  evident.  (C) 

Conclusion 


(C)  44.  Analysis  of  data  from  small-  and  large-yield  explosion  experi¬ 
ments  indicates  that: 

a.  For  depths  of  burial  greater  than  0.20  ft/lb  '  ^ 

(~25  ftAt1/3)  HE  and  NE  are  approximately  of  equal  effec¬ 
tiveness;  for  shallower  burial,  an  equivalence  factor  must 
be  used. 

b.  Large-scale  surface  explosions  can  be  modeled  using  confined 
charges  and  an  equivalence  factor. 
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_c.  Model  experiments  can  be  conducted  by  either  scaling  mate¬ 
rial  properties  or  using  the  same  material  in  the  model 
as  exists  in  the  prototype. 

d.  Lithostatic  pressure,  which  can  be  modeled  artificially, 
produces  significant  stress  for  depths  greater  than  3000  ft; 
therefore,  the  lithostatic  stresses  should  not  be  neglected 
when  attempting  to  model  explosion  effects  on  structural 
inclusions  that  are  deeper  than  3000  ft  in  the  prototype. 

e.  Scaled  closure  distances  agree  closely  with  the  cube-root 
law  with  some  allowance  made  for  variation  in  material 
properties. 

f.  Peak  strain  multiplied  by  pulse  duration  scales  as  the  cube 
root  of  yield  for  given  scaled  ranges.  This  parameter  plus 
energy  density  are  believed  to  be  the  important  ones  for 
design  purposes.  (C) 

(C)  45.  The  foregoing  analysis  indicates  that  model  scaling  can  be  used 
to  advantage  to  investigate  a  number  of  the  factors  which  are  important  in 
the  response  of  deep  underground  structures  to  explosive  attack.  (C) 

Recommendations 


(C)  46.  It  is  recommended  that: 

a.  The  present  investigation  be  pursued  by  experimentation 
both  with  models  which  scale  the  properties  of  the  material 
and  those  which  utilize  material  in  the  model  which  has 
the  same  properties  as  the  prototype. 

b.  Model  studies  of  the  effect  of  geologic  structures  be 
investigated. 

_c.  Attempts  be  made  to  model  lined  tunnels  w?th  the  linings 
backed  by  various  packings,  as  well  as  to  use  models  to 
study  vibration  effects  on  spring -mounted  structures. 

d.  Factors  related  to  modeling,  such  as  NE-HE  equivalence, 
effects  of  shallow  burial,  and  coupling,  be  investigated 
to  permit  more  extensive  extrapolation  of  model  results,  (c) 
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Adapted  from  reference  8  of  Selected  Bibliography. 

With  venting  to  the  unusable  portion  of  the  drift  beyond  the  zero  point. 
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Table  2 

Model  Ratios  of  Mechanical  Quantities* 


Scaling  Parameters 

(1) 

Dimensional 

Formula 

(2) 

Model 

Patio 

(3)** 
a  *  a 
>  P  a 

pP  =  p2m 
xVt>  =  V2m 

Acceleration 

LT~2 

vr'2 

1 

X-1 

Action 

ML2T_1 

pA2T-1 

mx3/2 

k 

X 

Angle 

L° 

1 

1 

1 

Angular  acceleration 

t“2 

-2 

T 

x"1 

X-2 

Angular  momentum 

ml2t_1 

HX2t_1 

MX^2 

X4 

Angular  velocity 

T_1 

T-1 

X"1/2 

X-1 

Area 

L2 

x2 

x2 

x2 

Curvature 

L'1 

X-1 

X-1 

x"1 

Density 

ML"3 
-1  -2 

MX’3 
-1  -2 

MX"3 

-2 

1 

Elastic  modulus 

ML  T 

-2 

|iX  T 
-2 

MX  * 

1 

2 

Force 

MLT 

|iXT 

M 

X 

Frequency 

t"1 

-1 

T 

x'1>/2 

x-1 

Gravitational  constant 

m_1l3t"2 

-13  -2 

M  XT 

M-1X2 

x-2 

Impedance  (pc) 

ml'V1 

-2  -1 
(iX  T 

MX'/72 

1 

Kinematic  viscosity 

lV1 

,2  -1 

X  T 

x3/2 

X 

Length 

L 

X 

X 

X 

Moment  of  momentum 

2  -1 

ML  T 

^2 
|iX  T 

MX372 

x1* 

Momentum 

MLT”1 

x  -1 

|iXt 

MX172 

X3 

Power 

ml2t"3 

^2  -3 

pX  T 

MX172 

x2 

Strain 

L° 

_ -1-2 

1 

-1  -2 

1 

-? 

1 

Stress ,  pressure,  and  strength 

ML  I 

MX  t 

ro 

1 

Time 

T 

T 

X 

Torque 

2  -2 

ML  T 

,2  -2 

MX  T 

MX 

X3 

Total  impulse 

MLT'1 

.  -1 
pXT 

.1/2 

MX  7 

X372 

Total  impulse  area 

ML-1T-1 

mx'V1 

MX’372 

X172 

X 

Velocity 

LT"1 

XT'1 

1 

Viscosity 

ML-1!'1 

£ 

l 

H 

MX"372 

X 

Volume 

L3 

X3 

X3 

X3 

Wavelength 

L 

X 

X 

X 

Work  and  energy 

2  -2 

ML  T 

.2  -2 
pX  t 

MX 

X3 

*  See  glossary  for  definition  of  symbols. 

**  a  ,  p  ,  and  V  stand  for  acceleration,  density,  and  velocity,  respectively 
The  subscripts  p  and  m  stand  for  prototype  and  model,  respectively. 
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Crater  Radius  Data  From  Full-Scale  HE  and  ME  Shots 
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SEDAN  2  x  10$  585  640  633  1. 09  1. 08 

TEAPOT  ESS  2. 4  X  10$  134  146  67  1. 09  0. 50 

JANGLE  S  2. 4  x  10$  134  45  -3. 5  0. 34  -0. 026 

JANGLE  U  2.4x10$  134  129  17.0  O.96  0.127 

JOHNIE  BOY  1. 0  x  10°  100  6l  2.7+  0. 6l  0. 03 
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Table  4 

Apparent  Volume  and  Radius  Efficiencies  of 
Full-Scale  HE  and  ME  Shots 

charSe  ~  r/W1^3  v/Wl/3 

Weicht  c_  a  >  V" 


Project 

weignp 

lb 

ft/lb1^3 

ft/lb1/3 

JL 

ftVib1/3 

t 

SCOOTER  (HE) 

1 

X 

1# 

1. 25 

1. 54 

100 

2.64 

100 

SEDAN 

2 

X 

108 

1. 08 

i.n 

72 

0.90 

34 

TEAPOT  ESS 

2.4 

X 

106 

0.50 

1.09 

71 

1.08 

40.8 

JANGLE  S 

2.4 

X 

106 

-0.  026 

0.34 

22 

1.85  x  10"3 

0.07 

JANGLE  U 

2.4 

X 

106 

0.  0127 

0.96 

62 

0.42 

16 

JOHNIE  BOY 

1 

X 

106+ 

o.  03+ 

o.6l 

4o 

0.15 

6 

Table  5 

Properties  of  UET  and  HAEDHAT  Rocks 


Property 

Unaweep 

Granite 

Navajo 

Sandstone 

HAEDHAT 

Granite 

Tensile  strength,  psi 

500 

100 

10,000 

Compressive  strength,  psi 

25,000 

4800-8600 

25,000 

Young's  modulus  of  elasticity, 
psi 

3. 4-5. 4  x  106 

1. 5-2. 5  X  106 

11.  3  x  lo1 

Modulus  of  rupture,  psi 

2,300 

600 

1,900* 

Strain  at  compressive  fail¬ 
ure,  psi 

0.004-0.007 

0. 003-0. 006 

0. 002 

*  Splitting  strength. 
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